Stroboscopic model of transport through a quantum dot with spin-orbit scattering 
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We present an open version of the symplectic kicked rotator as a stroboscopic model of electrical 
conduction through an open ballistic quantum dot with spin-orbit scattering. We demonstrate 
numerically and analytically that the model reproduces the universal weak localization and weak 
anti- localization peak in the magnetoconductance, as predicted by random-matrix theory (RMT). 
We also study the transition from weak localization to weak anti-localization with increasing strength 
of the spin-orbit scattering, and find agreement with RMT. 
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I. INTRODUCTION 

Electrical conduction in semiconductor heterostruc- 
tures is affected by the spin degree of freedom through 
spin-orbit scattering. In quantum dots with chaotic scat- 
tering a statistical approach is appropriate. The struc- 
ture of the spin-orbit Hamiltonian (of either Rashba 
or Dresselhaus form) has a special structure, that of a 
non-Abelian vector potential. By a gauge transforma- 
tion Aleiner and Fal'ko identified all possible symmetry 
classes and described the crossovers between them by 
means of random-matrix theory (RMT). 1 This RMT has 
been extended by Brouwer et al. to the case that the spin- 
orbit scattering is nonuniform and thus the gauge trans- 
formation cannot be made. 2,3 Experiments are in good 
agreement with the predictions of the theory. 4,5 Recently 
a semiclassical theory of quantum dots with spin-orbit 
scattering has been developed. 6 ' 7 

Here we present a fully quantum mechanical computer 
simulation to test the theory. In the case of spinlcss 
chaotic quantum dots, the stroboscopic model known 
as the quantum kicked rotator has been proven to be 
quite successful. 8-14 This model exploits the fact that, 
although the phase space of the open quantum dot is 
four-dimensional, the dynamics can be described, on time 
scales greater than the time of flight across the dot, as a 
mapping between points on a two-dimensional Poincare 
surface of section. The kicked rotator gives a map on 
a two-dimensional phase space that has the same phe- 
nomenology as open quantum dots. 

In this paper we extend the model of the open kicked 
rotator to include spin-orbit scattering in a perpendic- 
ular magnetic field. We begin by describing the known 
model for a closed chaotic quantum dot 15 with spin-orbit 
scattering in Sec. II A, before discussing the opening up 
of the model in Sec. II B. The relation of the model to 
RMT is given in Sec. III. This relation will give us a map- 
ping between the model parameters and the microscopic 
parameters of a chaotic quantum dot. Numerical results 
for the weak (anti)-localization peak and its dependence 
on magnetic field and spin-orbit scattering strength are 
presented in Sec. IV and compared with the analytical 



predictions from Sec. III. 



II. DESCRIPTION OF THE MODEL 

A. Closed system 

The symplectic kicked rotator has been introduced by 
Scharf 15 and studied extensively in Refs. 16-18. We sum- 
marize this known model of the closed system before pro- 
ceeding to the open system in the next subsection. 

The symplectic kicked rotator describes an electron 
moving along a circle with moment of inertia Iq, kicked 
periodically at time intervals tq with a kicking strength 
that is a function of position and spin. We choose units 
such that tq = 1 and h = 1. The Hamiltonian H is given 
by 15 - 16 

H=Up + p ) 2 + V{9) jr 6 s (t-n), (2.1a) 



n— — oo 



V{6) = K cos(0 + O ) + K so (ai sin 29 + a 3 sin 9). 

(2.1b) 

We have introduced the symmetrized delta function 
S s (t) = [S(t + e) + S(t — e)]/2, with e an infinitesi- 
mal. The dimensionless angular momentum operator 
p = —i% e gd/d9, with h e g — %tq/Iq the effective Planck 
constant, is canonically conjugate to the angle 9 S [0, 2tt). 
The kicking potential V(9) contains the Pauli spin ma- 
trices 



o\ = 



1 

1 



(T 2 



-i 

1 



CT3 = 





-1 



,(2-2) 

Potential scattering is parameterized by the kicking 
strength K and spin-orbit scattering by K so . 

Spin rotation symmetry is broken if K so =/= 0. The 
generalized time-reversal symmetry 15 



T 



(2.3) 



is preserved if 9a = and is broken if 9q E (0, ir). A 
nonzero po ensures that the Hamiltonian has no other 
unitary or antiunitary symmetries. 15 
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Notice that the roles of p and 8 are interchanged in T 
compared to the conventional time-reversal symmetry of 
the Rashba Hamiltonian and the spinless kicked rotator, 
which reads 



T : 8^8, p> 



-P, <?i 



t ■ 



-t. 



(2.4) 



For this reason time-reversal symmetry in the symplectic 
kicked rotator is broken by a displacement of 8, rather 
than by a displacement of p as in the spinless kicked 
rotator. 19 

The stroboscopic time evolution of a wave function is 
governed by the Floquet operator 



T = T exp 



IcS Jt 



to+l 



H(t)dt 



(2.5) 



where T denotes time ordering of the exponential. In the 
range [—1/2,1/2) only to ~ and to = —1/2 preserve 
T-symmetry for do = 0. We will find it convenient to 
choose to = —1/2 for numerical calculations and t = 
for analytical work. 

The reduction of the Floquet operator to a discrete 
finite form is obtained for special values of h e s, known as 
resonances. 19 For h c s — 4n/M, with M an integer, the 
Floquet operator is represented by an M x M matrix of 
quaternions. (A quaternion is a 2 x 2 matrix q = qol + 
iqi<Ji + iq2<?2 + iq3&3 with a complex number and 1 the 
2x2 unit matrix.) For this value of h e g the momentum is 
restricted tope [0, 47r), i.e. one can think of the Floquet 
operator as describing a map on a torus. For t = —1/2 
the matrix elements in the p-representation are given by 



T w = (UUXU*IL)w, U' = 0,1,...,M-1, (2.6a) 
U u .=6 w e-*W' U 1, l P ° M 

U w = M -^ 2 e- i2nU '/ M t, 
X w = <5 w , e - i(M/47r)v(27ri/M) . 



47T 



(2.6b) 

(2.6c) 
(2.6d) 



For to = one has instead 

T =UX l ' 2 U^U 2 UX l ' 2 U^. 



(2.7) 



The generalized time-reversal symmetry (2.3) is ex- 
pressed by the identity 



T = T R , if do = 0. 



(2.8) 



The superscript R denotes the dual of a quaternionic 
matrix, 

T R = a 2 T T a 2 . 

Here T denotes the transpose in the basis of cigenstates 
of p (p- representation) . To verify Eq. (2.8) note that 
= and that the transpose in p- representation 
takes 8 to —6. 




FIG. 1: Location of the absorbing boundary conditions (grey 
rectangles) in the classical phase space of the open kicked 
rotator. To ensure that the openings do not break the time 
reversal symmetry they are oriented parallel to the p-axis in 
the spinless kicked rotator (left panel) and parallel to the 8- 
axis in the symplectic kicked rotator (right panel). 



B. Open System 

To describe electrical conduction we open up the kicked 
rotator, following the general scheme of Refs. 8-11. We 
model a pair of iV-mode ballistic point contacts that cou- 
ple the quantum dot to electron reservoirs, by impos- 
ing open boundary conditions in a subspace of Hilbcrt 
space represented by the indices The subscript 

k = 1,2, ...TV, with N = Ni + N 2 , labels the modes 
(both spin directions), and the superscript fi = 1, 2 labels 
the point contacts. The N x M quaternionic projection 
matrix P is given by 



Pkk' 



if k' 
otherwise 



(m) 

k ' 



(2.9) 



The matrices P and T together determine the scatter- 
ing matrix 



S(e) = P(e- ie - TQ T QY 1 TP T , 



(2.10) 



where e € [0, 27r) is the quasi-energy and Q T Q = 1 — 
P T P. One readily verifies that S is unitary. 

We need to ensure that the introduction of the point 
contacts does not break the T-symmetry 



S(e) = S R (e), if <9 = 0, 



(2.11) 



or for non-zero #o the more general duality relation 

S(8o) = S R (-8 ). (2.12) 

This is assured by choosing the absorbing boundary con- 
ditions in a strip parallel to the 0-axis, rather than par- 
allel to the p-axis as in the spinless kicked rotator (cf. 
Fig. 1). The difference is due to the exchange of the 
roles of p and 8 in the time-reversal symmetry operation, 
compare Eqs. (2.3) and (2.4). 

By grouping together the indices belonging to the 
same point contact, the N x N quaternionic matrix S can 



3 



be decomposed into 4 sub-blocks containing the quater- 
nionic transmission and reflection matrices, 



S = 



r t 
t' r> 



(2.13) 



The value of s is arbitrary; we will take e — in the 
analytical calculations and average over e in the numer- 
ics. The T-symmetry (2.11) requires that r = <j 2 r 7 <j 2 , 
r' = <J 2 r' T a 2 , and t' — a 2 t T (j 2 . 

The conductance G follows from the Landauer formula 



G = — Tr tt\ 
h 



(2.14) 



where the trace Tr is over channel indices as well as spin 
indices. Unitarity of S ensures that Tr ift = Tr t'i't. For 
#0 = the eigenvalues of ttf are doubly degenerate due 
to the T-symmetry (Kramers degeneracy). It will prove 
useful to write the Landauer formula in the form 2,3 



G 



2e 2 NiN 2 



— Tr5A5 t A = G Q + SG, (2.15) 



h N h 
with A a diagonal matrix having diagonal elements 



Aw = 



N 2 /N 



1, 



,N lt 



—Ni/N j = m + l,...,N. 



(2.16) 



The term Go = (2e 2 /h)N\N 2 /N is the classical conduc- 
tance and the term 5G, of order e 2 /h, is the quantum 
correction from the weak localization effect. 



III. RELATION TO RANDOM-MATRIX 
THEORY 

Random-matrix theory (RMT) gives universal predic- 
tions for the quantum correction SG in Eq. (2.15). We 
calculate this quantity for the symplectic kicked rotator 
and compare with RMT. This will give us the relation 
between the parameters of the stroboscopic model and 
the microscopic parameters of the quantum dot. 

The three universality classes of RMT are labeled by 
(3 = 1,2,4, with 20 



SG 



RMT 



(3-2 e 2 
2/3 h 



(3.1) 



In the absence of T-symmetry one has f3 = 2. In the 
presence of T-symmetry one has (3 = 1 (4) in the presence 
(absence) of spin rotation symmetry. We will investigate 
the three symmetry breaking transitions (3 = 1 — > 2, 
(3 = 1 — > 4, and (3 = 4 — > 2 in separate subsections. 



(3 



2 transition 



The (3 = 1 — > 2 transition takes place in the absence 
of spin-orbit scattering (K so = 0). This transition was 



studied in Ref. 12 for the case that the symmetry T' 
rather than T is broken. To fully characterize the model 
we need to reconsider this transition for the case of T- 
symmetry breaking. 

For small 9 , cos(# + 9 ) w cos 9 — 8 sin9 and the 
Floquet matrix (2.7) takes the form 

T(K so = 0, 9 -> 0) = e e ° w F e 9aW , (3.2a) 
W = UYU\ Y w = Sun(KM/8n) sm(2nl/M). (3.2b) 

Here To = T(K so = 0, #0 = 0) is unitary symmetric and 
W is real antisymmetric. The scattering matrix (2.10) 
(at e = 0) becomes 



5 , = T(l-T i?)" 1 T T', 
T' = e eoM/ P T , 



(3.3a) 
(3.3b) 
(3.3c) 
(3.3d) 



Substitution of S into Eq. (2.15) gives the conductance 
G. 

To make contact with RMT we assume that To is a 
random matrix from the circular orthogonal ensemble 
(COE), expand the expression for G in powers of To and 
average To over the COE. In the regime 1 -c A^ -c M, 
we can perform the average over the unitary group with 
the help of the diagrammatic technique of Ref. 21. Since 
Tr A = only the maximally crossed diagrams contribute 
to leading order in N. The result for the average quan- 
tum correction becomes 

(SG) = -^tr ^KT { T>KT«r M _l R]RT . (3.4) 

The factor of 2 comes from the spin degeneracy and the 
trace tr is over channel indices only. The two remaining 
traces are evaluated in the limit N, M — > 00 at fixed 
N/M. We find 



M _1 tr T t AT(TAT' t ) T 



NiN 2 N 



(3.5) 



N 2 M ' 

M _1 tr R f R T = 1 - N/M - $1 {KM /An) 2 (1 - N/M). 

(3.6) 

Substitution into Eq. (3.4) gives the average quantum 
correction 



(SG) = 



e 2 2N\N 2 



1 



9r = 



h N 2 l + (9 /9 c ) 2 ' 
4ttV]V 



(3.7a) 
(3.7b) 



KM 3 / 2 ' 

The RMT result has the same Lorentzian profile 20 ' 22 
e 2 2A^ 2 1 



SG 



RMT 



B r = C 



h N 2 1 + (B/B c y 
h /NLA\ 1/2 



eh 2 \ Tivf 



■ 



(3.8a) 
(3.8b) 



4 



with C a numerical constant of order unity, L the size 
of the dot, A its mean level spacing and vp the Fermi 
velocity. Comparison of Eqs. (3.7) and (3.8) allows us to 
identify 



<W$c — B/B c . 



(3.9) 



/3 = 1 



4 transition 



The f3 = 1 — > 4 transition is realized by turning on 
spin-orbit scattering (-Kso) in the absence of a magnetic 
field (0 O = 0). In this transition the quaternionic struc- 
ture of the Floquet matrix plays a role. The Floquet 
matrix (2.7) has the form 



F(K SO ,6 = 0) = e K °° A f e K °° A , 
A = U(a 1 Y 1 + a 3 Y 3 )U\ 
(Yi)u> - -*«/i(M/87r)sin(47r//M) ) 
{Y 3 )w = -*«'i(M/87r)sin(27rf/M). 



(3.10a) 
(3.10b) 
(3.10c) 
(3.10d) 



The matrix A is real antisymmetric and thus A* = —A, 
where * denotes quaternion complex conjugation. (The 
complex conjugate of a quaternion q is defined as q* = 
q$ 1 + iq* <7i + iq% a2 + iq 3 03 • ) The scattering matrix takes 
the same form (3.3a), but now with 



T = Pe 



K so A 



K aa A D T 



R 



-,K 30 A 



Q T Qe 



K.„A 



(3.11a) 
(3.11b) 
(3.11c) 



The average of T§ over the ensemble of unitary sym- 
metric matrices only involves the channel indices and not 
the spin indices. To keep the quaternions in the correct 
order we adopt the tensor product notation of Brouwer 
et al. 2 ' 3 The average of SG over gives, to leading order 
in N, 



(SG) 



-E 
h ^ 

[IV 



tr E® E' 



Ml ® 1 - tr .R® R* 



(3.12) 



where r = 1(g) er 2 , E = T^AT, and E' = T'AT't. The ten- 
sor product has a backward multiplication in the second 
argument, 



(a<g>&)(c®d) 



• db, 



(3.13) 



and the indices /1 and v are the spin indices. 

The two traces are calculated in the limit K so — > 0, 
TV, M -» 00 at fixed N/M, leading to 



M _1 tr E®E h 
M~hv R®R* 



10 1, 



(3.14a) 



iV 2 M 

(1 - N/M)(l - 4 J ftT 2 (M/87r) 2 )l ® 1 



+ 2X 2 (M/8^) 2 (1 - iV/M)(<ri (g> ai + a 3 ® a 3 ). 



(3.14b) 



After substitution into Eq. (3.12) there remains a matrix 
structure that can be inverted, resulting in 



(5G) = 



e 2 ^1^2 



h N 2 
a = K so /K, 



1 



1 



1 



2a 2 1 + 4a 2 
47rv / 2iV 



M 3 /2 



(3.15a) 
(3.15b) 



The RMT result has the same functional form 2 , with 
a = (2ttHN y-TsoA) 1 / 2 . Here r so is the spin-orbit scattering 



(3.16) 



time. Thus we identify 

K so /K c = (2tt^/t so A) 1 / 2 . 

C. (3 = 4 — > 2 transition 



In the presence of strong spin-orbit scattering (K so 3> 
K c ) the Floquet matrix takes for small 0o the same form 
as in Eq. (3.2a) for K so = 0, but now = T(K so » 



CD 

to 



(5 




0.125 



FIG. 2: Average quantum correction {SG) to the conductance 
as a function of the T-symmetry breaking parameter 60 . The 
data points are for the symplectic kicked rotator character- 
ized by K = 41, M = 500, N = 10, l = 0.2. The solid lines 
are the analytical predictions (3.7) and (3.17) in the absence 
and presence of spin-orbit scattering. The dotted lines are 
the solid lines with a vertical offset, to account for a differ- 
ence between the predicted and actual value of the classical 
conductance Go- 



■5 




FIG. 3: Average quantum correction (SG) to the conductance 
as a function of spin-orbit scattering strength K so at zero mag- 
netic field. Other parameters are the same as in Fig. 2. The 
solid line is the analytical prediction (3.15) for the crossover 
from weak localization to weak anti-localization. 



uniformly distributed quasi-energies and 40 randomly 
chosen lead positions. The quantum correction (SG) is 
obtained by subtracting the classical conductance Go- 
The numerical data is shown in Figs. 2 and 3. The 
magnetic field parameter 9 is given in units of 9 C from 
Eq. (3.7b) and the spin-orbit scattering strength param- 
eter K so is given in units of K c from Eq. (3.15b). The 
solid lines are the analytical predictions (3.7), (3.15), 
and (3.17) without any fitting parameter. 

The small difference between the data and the predic- 
tions can be attributed to an uncertainty in the value 
Go of the classical conductance. A small vertical off- 
set (corresponding to a change in Go of about 0.1%) can 
correct for this (dotted lines in Fig. 2). The strongly non- 
Lorentzian lineshape seen by Rahav and Brouwer 13 ' 14 in 
the spinless kicked rotator is not observed here. 



V. CONCLUSION 



K c , 9 = 0) is a unitary self-dual matrix rather than a 
unitary symmetric matrix. We can repeat exactly the 
same steps as we did for K so = but with !F a random 
matrix in the circular symplcctic ensemble (CSE). We 
then average T$ over the CSE. This leads to 



(SG) = 



e 2 N^N 2 



1 



h n 2 i + (e /9 c y 



(3.17) 



with 9 C as in Eq. (3.7b). The width of the Lorentzian 
is therefore the same in the /? = 1 — > 2 and /3 = 4 — > 2 
transitions, in agreement with RMT. 20 



IV. NUMERICAL RESULTS 



We have presented a numerically highly efficient model 
of transport through a chaotic ballistic quantum dot with 
spin-orbit scattering, extending the earlier work on the 
spinless kicked rotator. Through a simple assumption 
of a random Floquct matrix we have derived analytical 
predictions for the conductance of the model as a func- 
tion of spin-orbit scattering strength and magnetic field. 
The functional form of the conductance coincides with 
random-matrix theory and through this correspondence 
we obtain a mapping from microscopic parameters to 
model parameters. Numerical calculations are in good 
agreement with the analytical predictions. This model 
may provide a starting point for the study of the effect 
of a finite Ehrcnfcst time on weak anti-localization, which 
requires very large system sizes (cf. Rcfs. 13 and 14). 



The numerical technique we use is the same as has 
been used before for the spinless kicked rotator. 11,12 A 
combination of an iterative procedure for matrix inver- 
sion and the fast-Fourier-transform algorithm allows for 
an efficient calculation of the scattering matrix from the 
Floquct matrix. 

The average conductance (G) was calculated with the 
Landaucr formula (2.14) by averaging over 60 different 
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